Phase diagram and critical points of a double quantum dot 
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We apply a combination of numerical renormalization group (NRG) and renormalized perturba- 
tion theory (RPT) to a model of two quantum dots (impurities) described by two Anderson impurity 
models hybridized to their respective baths. The dots are coupled via a direct interaction Uw and 
an exchange interaction J. The model has two types of quantum critical points, one at J — J c 
to a local singlet state and one at U12 = U12 to a locally charge ordered state. The renormalized 
parameters which determine the low energy behavior are calculated from the NRG. The results 
confirm the values predicted from the RPT on the approach to the critical points, which can be 
expressed in terms of a single energy scale T* in all cases. This includes cases without particle-hole 
symmetry, and cases with asymmetry between the dots, where there is also a transition at J = J c . 
The results give a comprehensive quantitative picture of the behavior of the model in the low energy 
Fermi liquid regimes, and some of the conclusions regarding the emergence of a single energy scale 
may apply to a more general class of quantum critical points, such as those observed in some heavy 
fermion systems. 
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I. INTRODUCTION 

Experimental evidence in heavy fermion materials sug- 
gests that there is a close connection between the regions 
of anomalous spin fluctuations associated with a quan- 
tum critical point and the occurrence of superconduc- 
tivity. There is a similar association in the cuprates 1 
and iron pnictide 2 superconductors. This has moti- 
vated many experimental studies of quantum behavior 
in fermionic systems, which can be induced by lowering 
the transition in a magnetically ordered state to zero by 
application of pressure, a magnetic field or by alloying^—. 
Many conjectures have been put forward to understand 
this behavior, but as yet there is no fully quantitative or 
comprehensive explanation (for a recent review on this 
topic see reference^). This situation has led to a resur- 
gence of interest in local or impurity models which have 
quantum critical points. There are two reasons for this. 
Firstly, impurity problems are simpler and well estab- 
lished techniques have been developed for tackling them. 
Secondly, such systems can be simulated using nanoscale 
devices, such as quantum dots, where the parameters can 
be varied in a controlled way via gate voltages and the 
behavior examined over different regimes. An example 
is the two channel Kondo model, which has a quantum 
critical point and non-Fermi liquid behavior. An arrange- 
ment of quantum dots and gate voltages has been used 
successfully to test experimentally theoretical predictions 
for this model in the non-Fermi liquid regime 5 . 

The two impurity Kondo model is another example of 
a model with local quantum critical behavior. Numer- 
ical renormalization group studies for this model have 
shown that there is a discontinuous transition on increas- 
ing the (antiferromagnetic) exchange interaction between 
the two impurities from a state in which the impurity 
spins are predominantly screened by the conduction elec- 
trons (Kondo screening) to one screened within a local 
singlet stated. There has been recent experimental work 8 



to examine the critical behavior of this model using two 
magnetic impurities, one on a STM tip and the other on a 
metal surface, such that the interaction between the two 
impurities can be modified by controlling the distance of 
the STM tip from the surface. Yet another type of criti- 
cal point which has been studied theoretically^, one that 
could occur in two capacitatively coupled quantum dots, 
has a discontinuous transition to a locally charge ordered 
state. 

In a recent paper we studied a two-impurity Ander- 
son model which has direct and exchange interactions 
between the two impurities and displays both types of 
transition^. We showed, by using a combination of 
renormalized perturbation theory (RPT) and numerical 
renormalization group calculations (NRG) , that we could 
derive many exact results for the low temperature behav- 
ior in the Fermi liquid regime right up to the quantum 
critical points. We were able to show the emergence of 
a single energy scale on the approach to each type of 
quantum critical point with renormalized Fermi liquid 
parameters corresponding to a strong correlation regime. 
The results were restricted to a model with symmetry 
between the two impurities and particle-hole symmetry. 
In this paper we enlarge upon those results, examine the 
phase diagram more fully, and also relax the restrictions 
to include channel asymmetry and lack of particle-hole 
symmetry. 
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FIG. 1. (Color online) A schematic illustration of NRG chains 
for the double dot system. The quantum dots (filled squares) 
interact via an exchange J and a direct term U12, and are 
hybridized with the first site of their respective conduction 
electron chains (filled circles). 
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The model for the system we study corresponds to two 
quantum dots connected via leads with their respective 
conduction baths. The Hamiltonian H for the system 
takes the form T~L = Hi + TL2 + TL12, where the individual 
dots are described by a single impurity Anderson model, 
so H a for dot a = 1, 2 is given by 



A'.IT 



where cft a a , d a a , are creation and annihilation opera- 
tors for an electron at the impurity site in channel a, 
where a = 1,2, and spin component a =t, i- The 
creation and annihilation operators, c\ a a , c k , are 
for conduction electrons with energy t k a in channel 
a. The hybridization function for each dot is given by 
A a ^{uj) = 7r^ fe \Vk, a \ 2 5{uj — efc). We will make the usual 
approximation and take A Q a (u>) to be a constant A Q , in- 
dependent of lo and cr, corresponding to the wide band 
limit with a constant density of states. 

The interaction between the two quantum dots is taken 
in the form of a direct term U12, and an exchange term 
J, so that the interaction Hamiltonian H12 nas the form, 



U12 = U 



12 



+ 2JSd.i • SrfO) 



(2) 



where J > for an antiferromagnetic coupling and J < 
in the ferromagnetic case. 

Though the model here is taken to describe a dou- 
ble quantum dot, it can serve equally well to describe a 
number of other physical situations. It can describe a 
single magnetic impurity with two fold 'orbital' degener- 
acy, where the exchange coupling J then corresponds to 
a Hund's rule term where a ferromagnetic exchange term 
would be appropriate^—. In the strong coupling regime 
U/ttA 3> 1, the model corresponds to two Kondo models 
and with U12 = it has been much studied as a model 
of a system with two coupled magnetic impurities^^— . 
The model with J = has also been used to study two 
capacitatively coupled quantum dots 9 . 

For the numerical renormalization group calculations 
(NRG) the conduction electron states are transformed 
to a basis corresponding to a tight-binding chain such 
that the model takes the form illustrated in Fig [TJ The 
densities of states for the conduction electron baths (half- 
bandwidth D) are discretized with a parameter A(> 1) 
as in the original calculation of Wilson^ so that the cou- 
plings along the chain fall off as A~ N / 2 . for large N, where 
N is the Ath site along the chain from the dot. 



II. PHASE DIAGRAM 

As the model used here is rather more general than 
those studied so far, we show some results for the phase 



diagram which now includes two QCPs. Before consider- 
ing the general case we look at the case with U12 = 0. In 
the Kondo regime this is a much studied model of a sys- 
tem with two coupled magnetic impuritie a ' 14 ~ ^£ which 
has a QCP at J = J c where there is a sudden change 
in phase shift of the conduction electrons from 7r/2 to 
zero. When J = the two dots are isolated and their 
ground states are (Kondo) singlets giving a phase shift 
of 7r/2 in each channel. In the opposite limit J — > 00 the 
impurity spins form a decoupled local singlet such that 
the phase shift of the conduction electrons in each chan- 
nel is zero. In some related models the change between 
these two limits was found to occur continuously with 
increase of J. Affleck, Ludwig and Jone o 18 i 19 , however, 
cleared up the confusion resulting from apparently con- 
flicting results from the different models, in giving the 
conditions for the change to occur discontinuously, re- 
sulting in a QCP. In a more recent study of a quantum 
dot version of the model, it has been shown that the tran- 
sition occurs even when there is asymmetry between the 
two channels and away from particle-hole symmetry 20 . 
For the calculations of the phase diagram we assume 
both channel symmetry Ui,V\ = U2, V2 and particle- hole 
symmetry Cd, a = —U a /2, though we relax some of these 
conditions later. It will be useful to introduce a renor- 
malised quantity A(0), which will be defined more fully 
later, such that A(U, 0) — > A for U -» and in the 
Kondo regime U/nA > 1, ttA(J7,0) -> 4T K , where T K 
the Kondo temperature of the isolated dot, defined by 
the relation Tk = (3/^B) 2 /4ximp> where Ximp is 

the T = 

spin susceptibility for an isolated dot. 
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FIG. 2. (Color online) A plot of J c /itA(U, 0) as a function of 
U/ttA for ttA = 0.01. 



In Fig. [2]we give a plot of the critical value J c /irA(U, 0) 
as a function of U/ttA for it A = 0.01 (note: we take this 
value for 7rA all subsequent calculations unless specified 
explicitly) . The transition in the NRG calculations corre- 
sponds to a sudden change in the low energy fixed point 
from that of a free chain with N sites, where N is even 
(odd) to one with N odd (even) . It can be seen that there 
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FIG. 3. (Color online) A plot of Ui 2 /ttA, as a function of 
U/nA for the model with J = 0. 



is a transition for all values of U including U = 0. For 
U/nA > 3 the curve is flat indicating that in this regime 
J c is proportional to Tk, corresponding to J c = 1.378Tk. 
These results are in line with calculations based on the 
Kondo model where values of J c /Tk ~ 1.4 were estimated 
when expressed using our definition of T ^ 14 ' 15 (in pre- 
vious studies Tk was defined by \s = 2/i|/7rTk giving in 
our case J c /Tk = 2.165). 

The second type of transition to a local charge or- 
dered state occurs in the model J = when U12 reaches 
a critical value f/f 2 - For U12 — the charges on the 
two impurities are identical, but on increasing U12 such 
that U12 3> U it becomes energetically unfavourable and 
there is a local broken symmetry resulting in a charge 
imbalance between the impurities. If there are matrix 
elements connecting the two possible broken symmetry 
states to restore the symmetry then they must be ex- 
ceedingly small because the NRG results indicate a pre- 
cise transition at a critical value U\2 = Vyi- We plot the 
value of U^/irA as a function of U/irA in Fig. [3] This 
transition exists for all values of U and for U /irA ^> 3 
the critical values approach the value of U, i.e. {7f 2 — ► U. 
These results are in agreement with those given earlier 
by Galpin et ali£. 

The phase diagram for the model with finite J and U12, 
where both types of transitions occur, is shown in Fig|3] 
for the value U/ttA — 5. The two types of transition 
are separated by the line U12 = U + 3 J/2 along which 
the phase shift remains at the value 6 = tt/2. The two 
types of transition asymptotically approach this line on 
opposite sides. 

To investigate the low energy behavior of the model in 
detail, we first derive a number of exact results in terms 
of renormalized parameters using a renormalized pertur- 
bation theory (RPT) for the low energy properties in the 
Fermi liquid regime. An analysis of the low energy NRG 
fixed point is then used to calculate the renormalized pa- 
rameters in terms of the 'bare' parameters that define 
the model. Once the renormalized parameters have been 



calculated, they can be substituted into the equations for 
the exact evaluation of the low energy thermodynamics. 
The relevant equations based on the renormalized pertu- 
bation theory (RPT) will be derived in the next section. 
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FIG. 4. (Color online) A plot of the critical lines (local singlet- 
stars, local charge order - circles) separating the regions with 
phase shift 8 = tt/2 and 8 = as a function of J/D and 
U12/D for U/nA = 5. The dotted line corresponds to U12 = 
U + 3J/2. The two vertical lines relate to results shown later 
in Fig. [26] and Fig. [27| 



III. RPT ANALYSIS 

We start with the Fourier transform of the single par- 
ticle Green's function for the impurity d-state, 

Gd,a,a{un') = - f <T T d a _ (r)d+ (0))e^ dr, (3) 



where w n > — (2n' + l)//3 and j3 = 1/T and the brackets 
(...) denote a thermal average. This Green's function can 
be expressed in the form, 



1 



iuj n - e d . a + iA Q>0 .sgn(a; n ) - S Q , ff (w n ) ' 

(4) 

where Y, a <T (uj n ) is the self-energy. 

For the zero temperature Green's function, which will 
be our main concern, u) n can be replaced by a continuous 
variable u>, and summations over u> n replaced by integra- 
tions over uj. For the renormalized perturbation theory^, 
the Green's function in Eq.(0| can be re-expressed as 

Gd, a ,a(oJ n ) = z a G d ^ a {u n ), where Gd, a , a (u) n ) is the 
quasiparticle Green's function given by 



Gd,a,a(^n) 



1 



iw„ - e d a + iA Q sgn(w n ) - S Q ^(w„) 



(5) 
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and the renormalized parameters, id a and A a are given 
by 



(6) 



where z a = 1/(1 — 3E Q , CT (cj)/9(iw)) evaluated at u = 0. 
The quasiparticle self-energy E Qi(T (w) is given by 



Sa,<r(w) 



2a S a , ff (u) - S O)O .(0) - 



<9£ Qj(T (w) 



diuj 



where we have assumed the Luttinger theorem 22 , 
ImS(O) = 0, so that ImS ai(7 (w) ~ w 2 as w — > 0. When 
expressed in this form, the uj — part of the self-energy 
and its derivative have been absorbed into renormalizing 
the parameters td,a and A Q , so in setting up the per- 
turbation expansion any further renormalization of these 
terms must be excluded, or it will result in over-counting. 
In working with the fully renormalized quasiparticles it 
is appropriate to use the renormalized or effective inter- 
actions between the quasiparticles. In the single channel 
case, we identify the renormalized interaction U a as the 
local four vertex r"j__j_ ^-(wi, 012,^3,^4) in the zero fre- 
quency Iimrfeai, 



(7) 



In this case we need to include the local inter-channel four 
vertex, T aa (cox , corresponds to the Fourier coefficient 
of the connected skeleton diagram for the two particle 
Green's function, 



(T T Sd,i(n) • S rf>2 (r 2 )), 
with the external legs removed. We define 

2i = 2 1 * 2 r oli0a (o,o,o,o). 



(8) 



(9) 



and we can define a renormalized direct interaction U12 
in a similar way. 

We can combine these terms to define a quasiparticle 
Hamiltonian TL = ^2 a T-L a + %12, where 



. e d,a 



di „d. 



+E 

k,a 

and 
T~Li2 



a^k,a,a 



ka 



U12 ■ 



nd,2.a' 



+2J : Sd,i ■ Sd,2 



( 10 ) 

The brackets .0: indicate that the operator O within 
the brackets must be normal ordered with respect to the 
ground state of the interacting system, which plays the 
role of the vacuum. This is because the interaction terms 
only come into play when more than one quasiparticle is 
created from the vacuum. 



The renormalized Hamiltonian is not equivalent to the 
original model, and the relation between the original and 
renormalized model is best expressed in the Lagrangian 
formulation, where frequency enters explicitly 23 . If the 
Lagrangian density C(ed, A a , U, J) describes the original 
model, then by suitably re-arranging the terms we can 
write 

C(e d , A, U, J, U 12 ) = C(e d , A a , U a , J, U l2 ) + C c ({Xi}), 

(11) 

where the remainder part C c ({Xi}) is known as the 
counter term. The set of coefficients {Ai}, where i = 
1,2,3,4,5, can be expressed explicitly in terms of the 
self-energy terms and vertices at zero frequency. These 
relations, however, are not useful in carrying out the ex- 
pansion. We want to work entirely with the renormalized 
parameters and carry out the expansion in powers of U a , 
J and Ui2- We assume that the Xi can be expressed 
in powers of quasiparticle interaction terms, U a , J and 
U12, and can be determined order by order from the con- 
ditions that there should be no further renormalization 
of quantities which have already been taken to be fully 
renormalized. These conditions are 



S a>CT (0) = 0, 



diuj 



0. 



(12) 



and that the renormalized local and intersite 4-vertices 
at zero frequency are unchanged. 

The propagator in the RPT is the free quasiparticle 
Green's function, 



4°Lk) = 



1 



e d>a + iA a sgn(uj r , 



(13) 



From Fermi liquid theory, the quasiparticle interaction 
terms do not contribute to the linear specific heat coeffi- 
cient 7 of the electrons. The specific heat coefficient 7 is 
given by 



7 = 27r 2 ]Tp Q (0)/3, 



(14) 



where p Q (w) is the free quasiparticle density of states per 
single spin and channel, 



A a /ir 



(15) 



{uj - i d , a ) 2 + A 2 a 
The results for the spin susceptibility is given by 

Xs = 2 lAz2'ns,aPa(0), Xc = 2^2vc,aP~ a (0), (16) 

a a 

where 

fjs.a = 1 + U a p~ a (0) - Jpp{0), 

Vc, a = l-U a p ct (0)-2U 12 pp(0), (17) 

where f3 ^ a. The equations for the spin and charge sus- 
ceptibilities, Xs and Xci can be proved using the Ward 
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identities which follow from conservation of total spin 
and charge. They also correspond to a simple mean field 
approximation to first order in the effective interactions, 
U a , J and U12, in the quasiparticle Hamiltonian Eqn. 
(TTUJ) . The corresponding equations for the staggered sus- 
ceptibilities also follow from these mean field equations 
and in the Appendix we give an alternative way they can 
be derived from the leading correction terms to the NRG 
fixed point, 



x:, 



where 



2^ 



2 V* ~st 



O«(o), 



xf =2£C**(°)' ( 18 ) 



vi% = i + u a p a (o) + JMQ)> 



(19) 



From the equation of motion of the Green's func- 
tion G^ ak , a {uS), which is the Fourier transform of 
the retarded conduction electron Green's function << 
c a k o-(^) • c a k' <r(*') ^ ^ m channel we can derive the 
equation, 



(w-efc) ' (w - e k ) " a ' a v (w-e fc /)' 

(20) 

where the t-matrix, i*'* (w) is given by t~'* (w) = 
14,aGd,Q !( T( w )14',a- From the definition of the phase shift 
as the phase of the t-matrix for u) = 0, 



tan 



_! / ImG d , a , g (0) 
VReG d , a , ff (0) 



(21) 



we find 



= —tan 



= —tan 



(22) 

The total occupation of the impurity sites ^ n.d jQ = 
5 a /^ at T = 0, which corresponds to the Friedel sum 
rule. 

We consider first of all the model with U\2 = J = 0, 
which corresponds to two independent Anderson models. 
With particle-hole symmetry €d. a = —U a /2, ed,a = 0, 
which corresponds to a phase shift 5 a — ir/2 in a given 
spin channel, independent of the value of U a . This leaves 
just two renormalized parameters per dot, A a and U a , 
which are related in the strong coupling regime, as there 
is only one relevant low energy scale, the Kondo temper- 
ature Tk,o- The relation takes the form, U a — nA a — 
4Tk jQ; , where Tk, Q is defined such that spin susceptibility 
of a single impurity is given by Xs.a = (3/iB) 2 /4T K , Q . 

On the approach to a quantum critical point we expect 
both quasiparticle weight factors, Z\ — > and z 2 — > 0, as 
a singularity develops in the self-energy of the dot Green's 
functions, which implies a divergence in the specific heat 
coefficient 7. It also implies a divergence in the quasi- 
particle density of states at the Fermi level, p a (0) — > 00. 
As J —¥ J c we expect a divergence in the staggered spin 



susceptibility xf as it is the fluctuations in this channel 
that become critical. We expect the other susceptibili- 
ties, Xs, Xc and to remain finite, which leads us to 
require rj SM — > 0, rj c ^ a — > and v^ a — > as J — » J c . 
These conditions give sufficient equations for us to deter- 
mine all the renormalized parameters /5 a (0), U a , J and 
U12 in terms of a single energy scale. We find 

U a p a [Q)^l, Jp Q (0)^2, U 12 p a (0)^0. (23) 

For the local charge transition it is only the the stag- 
gered charge susceptibility Xc which should diverge as 
U12 —> U\2- For the other susceptibilities to remain fi- 
nite, ?y SiQ — ► 0, ?]c. a —> and ry** Q — > as U12 — > L/f 2 . 
This implies 

U a p a (0)^-1, J/5 Q (0)^0, f/i 2 Pa(0)^l. (24) 

We can define an energy scale T* as J — > J c on the ap- 
proach to the QCPs via l//5(0) = 4T* which for particle- 
hole symmetry becomes irA = AT*, such that the renor- 
malized interaction parameters can be expressed in terms 
of this single energy scale T*. 
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FIG. 5. (Color online) A plot of A/A, U/irA and J/ttA as a 
function of J/Jc, for (7 = 0. 



Other exact results can be obtained from the RPT. An 
exact expression for the renormalized self-energy for the 
model with particle- hole symmetry £(w) to order oj 2 and 
T 2 follows from an RPT calculation 10 taken to second 
order in interaction parameters U, J and J7i 2 , and is 
given by 



E(w,T) 



2 IA 



G4 



(t*) 



TFT 



(25) 



where I = (2U 2 + 3 J 2 + 4l> 1 2 2 )/(7rA) 2 , so I -> 14 as 
J — > J c and I — > 6 as U — > £/f 2 , and we have taken the 
channel symmetric case, U\ = U2, V\ & = V%k,- This is 
a generalization of the calculations given earlier in a pa- 
per on an impurity model with a Hund's rule ter m 11 ' 12 . 
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FIG. 6. (Color online) A plot of the ratios, U/ttA and, and 
J /nA as a function of J/J c , for U — 0. 



On substituting (f25[) into the quasiparticle Green's func- 
tion, we deduce lu 2 and T 2 contribution to the t-matrix, 
t a (u,T), 



~\6 



(26) 



The ratio of the w 2 to the 7r 2 T 2 term is (I + 4)//, which 
for I = 14 gives 9/7 in agreement with the calculation of 
Mitchell and Sela 2 ^ for the two impurity Kondo model. 

To evaluate the RPT expressions for the low energy 
behavior of the model we need the values of the renor- 
malized parameters in terms of the 'bare' parameters that 
define the model. We give a brief description how they 
may be deduced from the low energy fixed point in an 
NRG calculation. The results of these calculations enable 
us to test the predictions given in Eqns (|23[) and (|2~4"|) . An 
evaluation of the RPT expressions gives a comprehensive 
picture of the low energy behavior of the model in the 
Fermi liquid regimes. 



FIG. 7. (Color online) A plot of Ai(U, J)/Ai(l7, 0) (dotted 
line and triangles), A2 (U, J)/ A2 (U, 0) (full line and circles), 
Ui/nAi(U,0) (dashed line and squares), U2/nA2(U, 0) (dot- 
dashed line and diamonds), and J/ir(Ai (U, 0) A 2 ((7, 0)) 1/2 
(full line and stars) as a function of J/J c , for [//VA2 = 6, 
A1/A2 = 4, A 2 = 0.01. 



3 

2.75 
2.5 
2.25 

2 

1.75 
1.5 
1.25 
1 

0.75 
0.5 
0.25 



1 








~~*--_ 










" 






r.H 






$ — $f — 

: 0— 


* 

-■e 


— » 


9 0.92 0.94 


0.96 


0.98 



1/1 



FIG. 8. (Color online) A plot of J/VAi (dots and 
squares), J/7rA 2 (dashes, stars), C/i/7rAi(full line, up trian- 
gle), Ui/n A2 (full line, down triangle), and A2/A1 (dot-dash, 
circles) as a function of J/J c for U/ti A2 = 6, A1/A2 = 4, 
ttA 2 = 0.01. 



IV. NRG CALCULATION OF RENORMALIZED 
PARAMETERS 

The numerical renormalization group results can be 
used to calculate the renormalized parameters, by iden- 
tifying the quasiparticle Hamiltonian in Eqn. (JTQJ) as 
the NRG low energy fixed point with the leading or- 
der correction terms. In the NRG calculation the non- 
interacting Green's function for the impurity site a takes 
the form, 



where 



(28) 



with n = 0, 1, 2, ....N for a tight binding iV-site conduc- 
tion chain, where V n are the intersite hopping matrix 
elements and e„ the energies, with V n = A~ n / 2 £„, where 
A (>1) is the discretization parameter, and is given 
by 



G 



(0) 

d,a,a 



1 



(27) 



D 



(1 + A- J )(l - A-"- 1 ) 



2 (1 - A- 2 »- 1 ) 1 /2(l _ A-2»-3)l/2- 



(29) 
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FIG. 9. (Color online) A plot of A/A, l d /nA, U/tyA and 
J/ttA as a function of J/J c , for U/tyA = 0.5, e d /^A = 0.159. 




FIG. 10. (Color online) A plot of l/p(0) (= 4T*) and e d p(0) 
as a function of J/J c , for U/ttA = 0.5, ed/vrA = 0.159. The 
discontinuity in e d p(0) at the quantum critical point corre- 
sponds to a jump in phase shift of ir/2 per spin state per dot 
channel. 



We use a discretization parameter A = 6, D = 1, particle- 
hole symmetry in the conduction band e n = 0, and retain 
4000 states at each NRG iteration. 

The rcnormalized parameters, id and V, are calculated 
by requiring the lowest energy single and hole particle ex- 
citations, Ep ta (N) and Eh, a (N) 7 for the interacting sys- 
tem with a NRG chain of length N to correspond to the 
poles of the quasiparticle Green's function in the limit of 
large N. This gives to the condition, 

E p/h , a (N) - e d - V 2 G ( ^ a {E p/Ka {N)) = 0, (30) 

,(0) 



where Gq U ^, a (uj) is given by Eq. This equation de- 

fines quantities, ld, a (N) and A a (N), which in general 
depend upon N. Plotting these quantities as a function 
of N el plateau develops when they become independent 
of N for large N, which then determines the parame- 
ters for the free quasiparticles. After diagonalizing the 



FIG. 11. (Color online) A plot of low energy local dot spectral 
density p(u)) = zp{ui) as a function of uj for the two cases 
J/Jc ~ 0.98 (full curve) and J / J c = 1-02 (dashed curve) for 
the same parameter set as in Fig. 1101 The inset shows p(0) 
(= z/5(0)) as a function of J/J c showing the discontinuous 
loss of spectral density at the Fermi level as the value of J 
increases through the critical point J = J c . 
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FIG. 12. (Color online) A plot of Up(0) and Jp(0) as a func- 
tion of J/J c , for U/tyA = 0.5, e d /TiA = 0.159. 



free quasiparticle Hamiltonian, the renormalized interac- 
tion parameters, U, J and U12, can be calculated from 
the leading order correction terms to the difference be- 
tween the lowest two particle excitation, E° p a (N) and 
the two lowest single particle excitations for large N, eg. 
E pp,a( N ) ~ E pA N ) ~ E p,ai N ) m channel a determines 
U a . Once the renormalized parameters have been deter- 
mined the T = susceptibilities and specific heat coeffi- 
cient 7 etc can then be determined by substituting them 
into the relevant RPT equations. For further details we 
refer to earlier paper a 12 ' 13 i 25 . This analysis differs from 
the one used by Krishnamurthy, Wilkins and Wilso n 26 i 27 
for the Anderson model where a distinction is made be- 
tween a strong coupling and weak coupling low energy 
fixed point. In our analysis there is just one type of low 
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energy fixed point for the Anderson model, the Fermi 
liquid one. The relation between the two approaches is 
discussed more fully in the Appendix. 

We now look in detail at the results for the renormal- 
ized parameters as the interactions between the dots are 
increased from zero to their critical values at the onset 
of the transitions. 
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FIG. 13. (Color online) A plot of A(U, J)/A(U, 0) as a func- 
tion of J/J c for values of U/ttA = 4, 5, 6. The results fall on 
a single curve indicating that this ratio becomes universal in 
the large U regime. 



A. NRG Results for U la = 0, J < J c 

We first of all show results of a calculation for the 
renormalized parameters for the symmetric model. Re- 
sults for U/irA — 5, which corresponds to the Kondo 
regime of the model, were given in earlier worki^, so here 
we show results in a contrasting regime with U = in 
Fig. \5\ Though U = 0, a renormalized interaction U is 
induced when the interaction J is switched on; though 
rather small initially, it increases until J ~ 0.7 J c and 
then falls to zero as J — > J c . It can be seen that all 
the renormalized parameters, 7rA, U and J tend to zero 
on the approach to the quantum critical point J = J c . 
The ratios, U /ttA and J/ttA, are shown in Fig. [51 The 
value of U/irA increases monotonically with increase of 
J and tends to the strong coupling value U / irA = 1 on 
the approach to the critical point. The ratio J/ttA also 
increases monotonically such that J/ttA — > 2 as J — >• J c . 
These results are in line with the RPT predictions in Eqn. 

We next look at a case where the two dots are hy- 
bridized differently to their respective baths so we have 
channel asymmetry. We take V\ — 2V 2 , which is such 
that Ai = 4A 2 , and take A 2 = 0.01 and U/nA 2 = 
6. For J = we are dealing with two isolated An- 
derson models where the ratio of the renormalized en- 
ergy scales of the two dots, Ai(J = 0)/A2(J = 0) = 



619.16. are very different. In Fig. [7] we give plots 
of A X {U, J)/A X {U,Q), A 2 (U, J)/A 2 (U,0), Ux/irAj&O), 
U 2 /ttA 2 (U,0), and J/7r(Ai(E7, 0)A 2 (U, 0)) 1 / 2 as a func- 
tion of J/J c . We see that U 2 = nA 2 (U,0) at J = 
showing that dot 2 is initially in the Kondo regime, but 
the ratio Ui/ttAx(U,0) < 1 so that initially dot 1 is not 
quite in the Kondo regime due to the larger value of the 
hybridization Ai (U, 0). As J is increased there is a tran- 
sition at J c /ttA 2 = 0.1605 where all the renormalized 
parameters tend to zero. The results imply that both 
z\ — > and z 2 — > 0, corresponding to a loss of the local 
quasiparticle weight at the Fermi level in both dots at 
the critical point. 

To test the predictions based on the RPT given in the 
Eqn. (|23p for this case, we plot the ratios of the renormal- 
ized parameters, J/nAi, J/ttA 2 , Ui/ttAi U 2 /ttA 2 and 
A2/A1 as a function of J/J c over the range 0.9 < J/ J c < 
1. Over this range we see that Ui/ttAi = U 2 /ttA 2 = 1, 
corresponding to strong coupling on both dots. Remark- 
ably the ratio A2/A1, which was initially very small, 
1/619.16 w 0.001615, approaches the strong coupling 
value 1 at the critical point. The two ratios, J/irAi 
and J/ttA 2 , converge to a common value 2 as J — > J c . 
These results confirm the predictions given in Eqn. (|2"5)) 
for the case U\ 2 = 0. Furthermore, they show rather 
dramatically the emergence of channel symmetry as J 
approaches the critical value J c . 

We found precisely similar results when we fixed the 
value of J at J/ttA 2 = 0.06 and U/irA 2 = 6, and in- 
creased the hybridization ratio r = A1/A2. A transition 
point was reached at r = r c = 2.6732, where all the 
renormalized parameters went to zero, and in the pre- 
dicted ratios. 

The transition in this model is remarkably robust and 
not only persists in cases where there is asymmetry be- 
tween the two dots, but also away from particle-hole 
symmetry 2 ^. Before looking at an example of this type 
we consider how to calculate the renormalized parame- 
ters for the regime J > J c . 



B. Extension of NRG results to J > J c 

1. Particle-hole asymmetric case 

For the regime J > J c there is a difference in the 
calculation of the renormalized parameters for the cases 
with and without particle-hole symmetry. In both cases 
the renormalized hybridization, which couples to dots to 
their respective conduction baths, V a — > as J — > J c for 
J < J c . In the particle-hole asymmetric case the renor- 
malized hybridization re-emerges in the regime J > J c , 
which means we can apply precisely the same analysis to 
the calculation of the renormalized parameters, A and 
£d, using Eqn. (I5H1) . as described earlier. To see what 
happens as we pass through the transition we look at the 
results for a particular example. 
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In Fig. [5] we show the results deduced for the renor- 
malized parameters over the range < J/J c < 1.2 
(J c /ttA = 1.0884) for the case with U/wA = 0.5 with 
tdj^A = 0.159, where the level on the dot is initially 
above the Fermi level. It can be seen that all the renor- 
malized parameters, A, id, U and J go to zero as J — > J c 
both for J < J c and J > J c . To learn in more detail what 
happens to the quantity e ( i through the quantum critical 
point we plot in Fig. [TUjthe product e<jp(0) and l/p(0) as 
a function of J / J c . There is a discontinuity in 6^/5(0) at 
J = J c which reflects the fact that there is a discontinu- 
ous change in the phase shift 6 a of ir/2. From the Friedel 
sum rule (|22[) this implies that the total local occupation 
number in the ground state for J > J c has two more elec- 
trons than that for J < J c . The local density of states 
p{u>) at each dot in the low frequency regime corresponds 
to zp(uj), and this is plotted in Fig. [TT]for cases on either 
side of the transition, J/J c — 0.98 and J/J c = 1.02. The 
peak in p(u>) has shifted from just above the Fermi level 
for J/J c = 0.98 to below for J/J c = 1.02. In the inset of 
Fig. [H]/9(0) (= z p(0)) is plotted as a function of J/J c . 
It can be seen that as the value of J increases through 
the transition there is a sudden loss of spectral weight at 
the Fermi level but there is finite spectral weight on both 
sides of the transition. 

To check the RPT predictions given in Eqn. (|2"51) we 
have plotted Up(0) and Jp(0) in Fig. [T2] It can be seen 
that Up(0) — > 1 and Jp(0) — > 2 as J — >• J c on both sides 
of the quantum critical point in agreement with the RPT 
prediction. 



2. Particle-hole symmetric case 

The case of particle-hole symmetry differs from the 
case we have just considered as the renormalized hy- 
bridization parameter V a , calculated from the NRG, does 
not just go to zero as J — > J c for J < J c but remains 
equal to zero for all J > J c . This implies that we lose 
all the spectral weight at the transition for the particle- 
hole symmetric case, and not just some of the spectral 
weight as was the case with particle-hole asymmetry. At 
the fixed point for J > J c the two dots are decoupled 
from the conduction electrons on the lowest energy scale. 
We know, however, for J > J c the low energy behav- 
ior still corresponds to a Fermi liquid. We conjecture 
that this Fermi liquid can be described by a similar set 
of renormalized parameters and that the equations for 
the susceptibility and low energy behavior are still valid. 
However, this needs to be clarified. 

The fact that z — > implies that the derivative of the 
the self-energy in the impurity Green's function, 
diverges at uj — such that this self-energy develops a 
singularity as J — > J c . The fact that z — for J > J c 
means that that this singularity persists in the regime 
J > J c and the analysis outlined in the previous section 
breaks down. To consider this situation in more detail, 
we sum over k and k' in Eqn. (|20p to obtain an equation 



for the Green's function on the first site of the conduction 
electron chain in the NRG calculation which we denote 
by G ,a :C rM, 



(31) 

If we take the conduction band in the wide band limit, 
V^Gq ^ a (oj) — —iirA a , this equation becomes 



(0) 



Go,«, ff (w) = (1 - iAG^M) • 



(32) 



If we introduce a self-energy £° a (uS) for the Green's func- 
tion for the first site on the a conduction electron chain, 
which will also include the effect of switching on the hy- 
bridization term, Eqn. (|32p becomes 



1 



i + ip c E£ i(r (w) = - — — — 

1 - lAG d ,a, 



(33) 



For J < J c the self-energies on the dot are non-singular 
and from the Friedel sum rule Ed + S Q (0) = corre- 
sponding to a phase shift S a = ir/2. We then have 
Gd, a .a{0) = 1/*A Q for the case of particle-hole symme- 
try and hence from Eqns. (|3"2")) and (|3"3")) the Green's 
function on the first site of the conduction chain is zero, 
Go,q.ct(0) = 0, and the self-energy E° a (uj) diverges at 
u = 0. 

For J > J c , however, the dot self-energy develops a 
divergence at to = 0, such that the dot Green's function 
vanishes for u> = 0, and we get a phase shift 5 a = 0. 
The self-energy for the first conduction site is now non- 
singular and from Eqn. (|3"3")) . E° ff (0) = 0. We can then 
develop a renormalized perturbation approach and anal- 
ysis of the low energy NRG fixed point with the first 
conduction site playing the role of an effective dot. We 
should, therefore, be able to describe the Fermi liquid 
behavior for J > J c in terms of renormalized parameters 
associated with the first site on the impurity conduction 
electron chain. 

For J < J c we determine the renormalized parameters 
id and V by requiring the lowest energy particle and hole 
excitations of the interacting systems to correspond to 
poles of the impurity quasiparticle Green's function. For 
J > J c we can adopt the same procedure but replacing 
the Green's function of the impurity site by that for the 
first conduction electron site, so that parameters for £q 
and Vo are required to correspond to solutions of the 
equation, 

E p/h , a (N) - e~ - V 2 G^ !(T (E p/hta (N)) = 0, (34) 

where Gi a a (u)) is given by Eq. (|2"5]l. The local quasi- 
particle interaction terms, U, J and U12, can then be 
calculated as in the case J < J c . 

As the effective impurities now correspond to the first 
site on the conduction electron chain for each channel, 
there are a some modifications to the analysis used for 
J < J c . The density of states of the conduction electron 
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chain which now starts at the second conduction site is 
modified, which affects the calculation of the effective hy- 
bridization width Ao = 7r|Vo| 2 Pc where p c is the density 
of conduction states at the Fermi level. The factor ^4a, 
which relates the hybridization width of the discretized 
model to the continuum model, also changes. These two 
effects can be combined into a single additional factor B\ 
so that the hybridization width Ao for the continuum 
model is calculated from A = 7r|Vb| 2 /2A^B^D. The 
value of B\ was estimated from the result for the isolated 
Anderson model (J = 0) for U/ttA > 1, U/ttA ->■ 1, 
which for a discretization parameter A = 6, gave the 
value B 6 = 2.2373. 

In the limit J — > oo the two impurities are entirely de- 
coupled from the conduction electrons. This is reflected 
in the NRG results for the renormalized parameter Vq, 
which in this limit give Vo — > A -1 / 2 ^, which is the value 
for the free NRG conduction chain. The calculation of the 
impurity susceptibility using the RPT expression with 
Ao includes this extra conduction electron contribution. 
However, this correction is negligible except in the regime 
when J / J c is very large and it is not necessary to take it 
into account. 




FIG. 14. (Color online)A plot of U/ivA and J/^A, as a func- 
tion of J/J c for U/ttA = 4, 5, 6. 



indicating universal behavior in this regime. We have 
dropped the index on Ao and Ao for the regime J > J c , 
to define a continuous function A through the transi- 
tion point. With the definition of the Kondo tempera- 
ture T K , nA(U, 0) = 4T K and the energy scale T*(U, J), 
irA(U, J) = 4T* we now have a single energy scale over 
the range < J/ J c < 2 such that in the strong correla- 
tion regime, 



T* =T K F(J/J C ), 



(35) 



where the function F(J/J C ) is independent of U for 
U/nA > 1, such that F(0) = 1 and F(l) = 0. We 
also have J c = 1. 3781k in this regime, so Eqn. (|35[) can 
be re-expressed in the form, T* = T K F (0.726 J/T K ). 

In Fig. Q3] we plot the ratios U/ttA and J/VA, for 
U/irA = 4,5,6 over the same range. All the values for 
U/irA fall on the same line U/irA — 1 corresponding to 
the strong coupling limit. The values of J/irA also all lie 
on a single curve. Hence, using Eqn. (|35[) , in this regime 
we have 



J = 4T*f(J/J c ), 



(36) 



where f(J/J c ) is a universal function such that /(0) = 
and /(l) = 2. As from Eqn. (l3"5j) T* is a universal 
function of J and Ik then Eqn. (|36|) can be re-expressed 
as J = Tk/(0.726J/Tk), where /(<// J c ) is a universal 
function such that /(0) = and /(l) = 0. 
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C. NRG results for J / 0, U w = 

We now consider the channel and particle-hole sym- 
metric model taking V± = Vi so we can drop the chan- 
nel index a. We extend the results for the renormal- 
ized parameters to include the regime J > J c . We 
restrict attention first of all to the strong correlation 
regime U/ttA 3> 3 for different values of U. In Fig. 
[13] we plot the results for the ratio A(U, J)/A(U, 0) as 
a function of J/J c for U/irA — 4,5,6 over the range 
< J/ J c < 2. It can be seen that, though the values of 
J c vary significantly, all three results fall on a single curve 



FIG. 15. (Color online) A plot of A(U, J)/A(U, 0) as a func- 
tion of J/J c for values of U/nA — 0,0.5, 1,6. 



The universality found in the Kondo regime U /it A ^> 
1 no longer holds for smaller values of U as can be 
seen in Fig. [15] where we give results for the ratio of 
A(U, J)/A(U,0) over the range < J/J c < 2 for values 
of U/ttA = 0,0.5,1,2,6 as a function of J/J c . It can be 
seen that they all qualitatively behave in a similar way, 
but for smaller values of U /irA the results do not lie on 
the same curve. 
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In Figs. [16] and [17] we give the corresponding plot 
for the ratios U/irA and J/nA. There is no universal- 
ity except in the approach to the critical point where 
U/irA — >■ 1 and J/irA — > 2 in all cases in line with the 
predictions in Eqn. 



FIG. 18. (Color online) A plot of the Wilson ratio fj B 
function of J/J c for U/irA = 4, 5, 6 for ttA = 0.01. 




FIG. 17. (Color online) A plot of J/nA, as a function of J/J c 
for U/wA = 0, 0.5, 1, 6. 



The factors fj s and fj c in the expressions for the spin 
and charge susceptibilities reflect the effect of the inter- 
action between the quasiparticles in enhancing or sup- 
pressing the corresponding fluctuations. For the single 
impurity Anderson model f) s is better known as the Wil- 
son ratio and in the strong coupling regime has the value 
2, enhanced over that for non- interacting quasiparticles 
where the value would be 1. In Fig. [18] we plot the values 
of fjs for U/irA = 4, 5, 6 as a function of J/J c . We see 
that there is a universal curve in this regime which falls 
from the value 2 for the isolated impurities to zero at the 
critical point. In Fig. [19] the values of fj s are plotted 
as a function of J/J c for U/irA = 0,0.5,1,4. There is 
non-universal behavior in the regime, though the results 



\ i 




FIG. 19. (Color online) A plot of fj s , as a function of J/J c for 
U/irA = 0, 0.5, 1, 4 for ttA = 0.01. 



The values of fj c for U/irA = 0, 0.5, 1 are shown in Fig. 
l20l For U = and J = the value is 1 corresponding to a 
non-interacting system and then, as J is increased, pro- 
gressively reduces to zero (corresponding to the Kondo 
limit) as J — > J c . For the larger values of U, fj c is al- 
ready reduced at J = 0, due to the on-site interaction U, 
and is then further reduced to zero as the exchange in- 
teraction increases to the critical value. For U/wA > 3.5, 
fj c ?s as the charge fluctuations are almost completely 
suppressed for all values of J in this range. 

All the expressions given in the previous section apply 
equally well to the situation with a ferromagnetic ex- 
change J < 0. In figure [21] we show the values of fj s as a 
function of J / nA extended to the ferromagnetic range for 
U/irA = 5, where it extrapolates to the expected value 
8/3. 
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[/i2/ttA, are shown as a function of U12/U over the same 
range. Over the range U12/U < 0.9 these ratios differ 
very little from the values for the Anderson model in the 
Kondo regime, U /irA = 1 and Ui 2 /irA — 0. Beyond this 
point there is a rapid change and the two curves cross at 
the SU(4) point where U/irA = U 12 /nA = 1/3. This re- 
sult follows from the RPT equations for the susceptibility 
from the condition that the uniform charge susceptibility 
\c is so small, due to the large value of U, that it can 
be equated to zero. Beyond the SU(4) point the ratios 
rapidly approach the values predicted in Eqn. (|24[) on 
the approach to the QCP, U /irA = — 1 and Ui 2 /ttA = 1. 
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FIG. 20. (Color online) A plot of f\ a , as a function of J/J c for 0.2 
U/ttA = 0, 0.5, 1 for ttA = 0.01. 
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FIG. 21. (Color online) A plot of fj s (Wilson ratio) as a 
function of J/nA, including the ferromagnetic range, for 
U/ivA = 5 for ttA = 0.01. 



D. 



Results for U12 < U? 2 , J = 



The transition to a locally charge ordered state in the 
model with a finite interaction U12 (J = 0) has been 
studied in detail by Galpin et al£. Here we throw further 
light on this transition by calculating the renormalized 
parameters. In particular we test the RPT conjecture of 
the emergence of a single energy scale on the approach 
to the QCP and the predictions given in Eqn. (f2"4")h 

In Fig. [55] we give results for the ratio A/A as a 
function of Un/U for the case U/ttA = 5. Over the 
range Un/U < 0.8 there is a slow steady increase in A 
with increasing U\2- There is then a rapid increase due 
to fluctuations of charge between the two dots and A 
reaches a maximum at the SU(4) point U12 = U, and 
an extremely rapid fall off beyond this point to a QCP 
at U12/U = 1.028. In Fig. [23] the corresponding re- 
sults for the renormalized parameter ratios, U/ttA and 
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FIG. 23. (Color online) A plot of U/nA (circles) and Uu/nA 
(stars) as a function of U12/U for U/nA — 5. The dotted 
line which passes through the SU(4) point corresponds to 
U 12 /nA = U/nA = 1/3. 
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The question arises as to what happens for smaller 
values of U and in particular the case U = 0. This is 
addressed in Fig. [24] which shows a plot of U/irA as a 
function of UyilTJy^ for values of U/irA = 0, 0.5, 1, 2, 3, 5. 
We see that in all cases U/irA — > — 1 as J7 12 — > U^, 
though less rapidly the smaller the value of U. This also 
holds for the case U = so that switching on the inter- 
action term U12 induces an onsite effective interaction U 
in the quasiparticle Hamiltonian. 

The corresponding curves for the ratio Ui 2 /irA as a 
function of U\2/U± 2 are shown in Fig. (23 We see that 
it this case all the curves approach the predicted value 
U^/irA — >• 1 as U\ 2 U12, though much less rapidly for 
the smaller values of U. 
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FIG. 24. (Color online) Plots of U /irA as a function of 
Uu/U^ for values of U/irA = 0, 0.5, 1, 2, 3,5, J = 0. 
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ues of both J and [/ 12 . We look first at a case with 
U/irA = U^/nA — 5 over a range of values of J through 
the critical point J c = 0.090157rA. This corresponds to 
moving along a vertical line in Fig. 0]with U12/D = 0.05. 
The results for the renormalized parameter ratios, J/irA, 
U/irA and U^/ttA are shown in Fig. [26l The results 
include the regime with ferromagnetic coupling J < 0, 
and it can be seen that as the ferromagnetic coupling 
increases U/irA approaches the value 1 and J/irA the 
value —2/3 as expected in the regime where the chan- 
nel and charge fluctuations are suppressed. The value of 
U^/irA is small in this limit but increases as the ferro- 
magnetic coupling is reduced. At J = there is an SU(4) 
point where the curves for U /irA and Uri/irA meet at a 
common value 1/3, and J = 0. As J is increased from 
J = the curves for J/irA and U/ttA increase to a maxi- 
mum at the transition point J = J c , where they reach the 
values 2 and 1, respectively. At this point U^/ttA = 
in line with the general prediction in Eqn. (I23[) . This 
means that switching on and increasing the value of J 
have the effect of cancelling out the [/ 12 /7rA ratio which 
is finite for J = 0, and reducing it to zero as J — > J c . For 
./ > J c , the values of J /n A, U/ttA fall off significantly 
with increase of J and U\2/nA has a slow increase. 

In Fig. [23 we plot the same ratios of the renormalized 
parameters but in this case we take value U\%/txA = 5.2. 
We take a range of values of J corresponding to moving 
along a vertical line in Fig. [4] U12/D — 0.052. We see 
in this case the line passes through three QCPs. For the 
ferromagnetic regime where (J — J c )/ J c < —6, the spin 
and inter-dot fluctuations are suppressed and the ratios 
take on the expected values, U/ttA — 1, J/irA = —2/3 
and Ui2/irA = 0. Around (J — J c )/J c ~ —4 there is a 
rapid change over as the values move to U/irA = — 1, 
J /it A — and U^/irA = 1 on the approach to the local 
charge order transition at (J — J c )/J c ~ —2.1. When 
(J — J c )/ 'J c ~ —0.25 the system emerges from the charge 
ordered state, with the same set of ratios, but as J — >• J c 
there is a rapid change of the ratios of the renormalized 
parameters to the ratios, U/irA = 1, J /it A = 2 and 
U12/1TA — on the approach to the transition to the local 
singlet state. There is a rapid decrease in U/irA, J/irA 
and an increase in the ratio U^/ttA. For (J — J c )/J c > 
0.25 there is a slow fall off with increase of J for all three 
ratios. 



FIG. 25. (Color online) Plots of U^/irA as a function of 
U12/U12 for values of U/irA = 0, 0.5, 1, 2, 3,5, J = 0. 



F. The case U12 = U + 3J/2 



E. NRG results for J ^ and U12 + 

To complete the check on the predictions in Eqns. 
and (f2~4"|) we look at some examples with finite val- 



There is a special line in the phase diagram shown 
in Fig. 21 corresponding to U\2 = U + 3 J/2 along which 
there appears to be no transition. With this set of param- 
eters, and a ferromagnetic exchange coupling (J < 0), the 
model corresponds to a model introduced by Yoshimori 11 
to describe a single impurity with an n-fold degenerate 
orbital for the case n = 2. The exchange interaction J 
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from Eq. ([37]) we find the value, |J| = 2ttA/3. The 
value fi s approaches rj s = 8/3, as expected for the Wilson 
ratio for a two channel Kondo model. 

We have previously calculated the renormalized pa- 
rameters for this model, both for the particle-hole sym- 
metric and antisymmetric models, for a ferromagnetic 
couplin g 12 ! 13 . It is interesting to check here what hap- 
pens in the antiferromagnetic case. The lack of a tran- 
sition is likely to be due to the fact that for the isolated 
dots at half filling the gain in energy in forming a local 
singlet —3 J/2 is compensated by an equal and opposite 
term from the J-dependence of the U12 term. We look 
at the results for the renormalized parameters in detail. 



FIG. 26. (Color online) A plot of the renormalized parameter 
ratios, J/nA, U/ttA and U\i/-kA as a function of (J — J c )/J c 
where J c = 0.09015ttA for U/tvA = U 12 /ttA = 5. As J is 
increased there is an SU(4) point at J = where U/irA = 
U\2/nA = 1/3 and J = followed by a transition at J = J c 
from the Fermi liquid phase J < J c to the second Fermi liquid 
phase J > J c . 
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FIG. 28. (Color online) A plot of the renormalized 
parameters, ttA/ttA, U /it A and J /it A as a function of J/nA 
for U/tvA = 5,U 12 = U + 3 J/2. 



FIG. 27. (Color online) A plot of the renormalized parameter 
ratios, J/nA, U/ttA and U\i/-kA as a function of (J — J c )/J c 
where J c = 0.27228ttA for U/tvA = 5, Ui 2 /ttA = 5.2. As J 
is increased there is first transition from a Fermi liquid phase 
J < J c to a local charged ordered regime indicated by the box 
(CO) followed by a return to the Fermi liquid phase J < J c 
and then a transition to the second Fermi liquid phase J > J c . 



in this context corresponds to a Hund's rule term. The 
model in the ferromagnetic regime was used to calculate 
the orbital susceptibilit y 12 ! 13 for a two-fold degenerate 
impurity with a Hund's rule interaction J, where the or- 
bital susceptibility is given by 



Xorb — 



i lorb = l + {U + 3J)p(0)- (37) 



As the value of \J\ is increased in the ferromagnetic range 
the orbital fluctuations are suppressed and the two local 
spins are coupled to give an effective spin 1 . When 7y or b ~ 
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FIG. 29. (Color online) A plot of ratios of the renormal- 
ized parameters, U /irA and J/ttA as a function of J/ttA for 
U/tvA = 5,Ui2 = U + 3J/2. 
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In Fig. [55] we plot the renormalized parameters 
ttA/ttA, U/ttA and J/ttA as a function of J/nA over an 
antiferromagnetic range of J for U/ttA = 5. The value 
of A, and consequently the quasiparticle weight factor 
z, decreases monotonically to very small values for large 
J without any evidence of a transition. Somewhat sur- 
prisingly in this case the value of U/ttA decreases and 
becomes negative with increase of J. 

In Fig. [55] we plot that ratios U/irA and J/ttA as a 
function of J/ttA. For J = 0, V /it A = 1/3, which is 
predicted from the condition that the f) c ~ due to the 
large value of U / ttA which suppresses the charge suscep- 
tibility. For large J, U/ttA -1/3 and J/nA -> 2/3, so 
there is a single energy scale for large values of J. These 
results are in line with predictions based on the fact that 
for large U and large J both the charge and spin suscep- 
tibilities are suppressed so fj c — > and t) s — ^ 0. From Eq. 
(|17p the condition fj c — >■ we find 3(17 + J) = ttA. Using 
this result in the expressions for fj s and fj or b we find 



2J 



??orb = o + 



2,7 

TTA' 



(38) 



From the condition /% — > 0, using Eq. (|17p . we find 
J/ttA — > 2/3. When this result is substituted into the 
earlier condition, 3(f7+ J) = irA, we find U/irA — ^ — 1/3 
in this limit, and fjorb — > 8/3. 

From earlier results for this model in the ferromag- 
netic ranged, when |J| was large such that the orbital 
susceptibility was suppressed we found J /it A — > —2/3 
and fj s — > 8/3, which suggests from Eq. (15"51) the rela- 
tion fj s o fjorb if J /it A o — J /it A under the sign change 
J — > — J. In Fig. [30] we plot the values of fj s , fj oro and 
fj c as a function of J/A(5, 0) for U/ttA — 5. The results 
show that there is an approximate symmetry fj s fj or b 
as J O — J, but it is not precise for values in the inter- 
mediate range of \J\ where the differences, though small, 
are greater than any that might arise from the estimated 
errors in the calculation of the renormalized parameters. 

In the antiferromagnetic regime there is still the ten- 
dency for the impurity spins to be screened more locally 
as J is increased. This can be seen in the plot of \s 
shown in Fig. [3T] The value of \s decreases monotoni- 
cally as J /it A increases to a value J/nA ~ 0.15. After 
that point there appears to be a slight upturn, but no 
real significance can be attached to this as it involves the 
product of a quantity tending to zero fj s and a diverging 
quantity 1/zirA as z — >• 0, and numerical errors start be- 
coming significant when z becomes very small. Though 
we get no transition in this case it does not correspond 
to the category described by Affleck et alJ^, where the 
phase shift S decreases with increases of J. For the whole 
range of J the phase shift remains at the value 6 = tt/2. 
It would appear that provided the Friedel sum rule ap- 
plies for particle-hole symmetry (rid^a ><r ) = 1/2 and if Eq. 
td + S(0) = 0, we always have a phase shift 6 = n/2, so 
we cannot have a continuous crossover from S = tt/2 to 
S = 0. Either there is no transition and S = tt/2 for all J 



or the self-energy develops a singularity at ui = 0, with a 
breakdown of the Friedel sum rule, and a sharp transition 
to a state with (5 = 0. 
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FIG. 30. (Color online) A plot of rj s (circles), f\ or b (squares) 
and 77 c (crosses) as a function of J/nA(5,0) for U12 — U + 
3J/2, U/ttA = 5. 




FIG. 31. (Color online) A plot of Xs(J)/Xs(0) as a function 
of J/ttA for U 12 = U + 3J/2, U/nA = 5. 



V. CONCLUSIONS 

We have examined the low energy behavior in the 
Fermi liquid regime of a model that has two types of 
quantum critical points. The thermodynamic and dy- 
namic response functions in this regime can be expressed 
in terms of a limited number of renormalized parameters, 
A, U, J and U12, which we can calculate accurately from 
an analysis of the NRG low energy fixed point. Once 
these have been determined they can then be substituted 
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into the relevant RPT formulae. The fact that certain 
susceptibilities remain finite at the transition, where the 
quasiparticle weight z — t 0, gives enough equations to 
predict the dimensionless parameters, Up(0), Jp{0) and 
f/i2/5(0) in the Fermi liquid region on the approach to 
the critical point. These predictions have been confirmed 
from the NRG results for both types of transition, includ- 
ing situations away from particle-hole symmetry, and dot 
asymmetry in the case of the transition at J = J c . As 
these dimensionless parameters are universal on the ap- 
proach to the critical points, the quasiparticle interac- 
tions, U, J and U12, can all be expressed in terms of a 
single energy scale T* = 1/4/5(0), where T* — > at the 
critical point. 

In this paper we have not examined the finite temper- 
ature non-Fermi liquid regime in the region of the crit- 
ical point at J = J c , but in our NRG results we found 
the same higher energy non-Fermi liquid fixed point as 
observed in previous calculations for the two impurity 
Kondo modet^. Those results were explained by Affleck 
and Ludwi g 18 ! 19 using conformal field theory. More re- 
cent work has clarified with the relation between this 
transition and the non-Fermi liquid fixed point of the 
two channel Kondo model. A recent paper of Sela and 
Mitchell^! has addressed the crossover as a function of 
temperature between the Fermi liquid and non- Fermi liq- 
uid regimes. It would be interesting to examine this 
crossover in terms of renormalized parameters using the 
approach which was developed and applied to the two 
channel Kondo model 28 -^. 

It has proved difficult to probe the quantum critical 
point of the two dot/impurity model experimentally. The 
problem is that in any experimental set up there will be 
in addition to the interdot/impurity interaction a direct 
or indirect hybridization term. This was the case in the 
recent experimental work on a two impurity system which 
consisted of a cobalt ion on a STM point and interacting 
with a second cobalt ion on a metal surface. The hy- 
bridization term destroys the critical point. Nevertheless 
it would be interesting to generalize the model used here 
to model this experimental system, to see if the quantum 
critical behavior could be observed due to the presence 
of a higher energy non- Fermi liquid fixed point, as was 
done in the case of the two channel Kondo model. 

The quantum critical points studied here are for impu- 
rity models and cannot be applied directly to the exper- 
imental results on quantum critical behavior observed in 
lattice heavy fermion systems. However, certain features 
found may be universal and apply to a general class quan- 
tum critical points. For example, the fact that on the 
approach to the quantum critical points all the quasipar- 
ticle interactions could be expressed in terms of a single 
energy scale T* means that any dynamic response func- 
tion should take the form xM = F(u/T*,T/T*). 
At the critical point where T* — > we would then ex- 
pect a form, xi 10 ) — T a f(ui/T), or u/T scaling which 
has been observed at the several heavy fermion critical 
points, such as in YbRu2Si2 30 . The quantum critical 



point induced in YbRu2Si2 by suppressing the antiferro- 
magnetic order with a magnetic field and has been inter- 
preted as a Kondo breakdown^, where the loss of the 
f-electron states at the Fermi surface results in a change 
in volume of the Fermi surface at the quantum critical 
point. A schematic sketch is given Fig. 4 in the paper of 
Pfau et alj^ of the evolution of the quasiparticle weight 
factor z across such a quantum critical point. At the crit- 
ical point J = J c in our calculations in the particle-hole 
symmetric case there is a similar Kondo breakdown be- 
cause the Kondo resonance at the Fermi level for J < J c 
collapses at J = J c and there are no impurity/dot f-states 
at the Fermi level for J > J c . In Fig. [13] we have pre- 
cise results for this evolution of the quasiparticle weight 
factor z for the Kondo collapse in a two impurity model 
(note that for comparison with the sketch for YbRu2Si2, 
the range J > J c corresponds to the small Fermi surface 
and J < J c to the large one). 
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VI. APPENDIX 

We clarify here how the analysis of the low energy NRG 
fixed points which we use in this paper differs from the 
approaches based on the work of Krishnamurthy, Wilkins 
and Wilson (KWW I and H) 26 i 27 . We consider the case 
of the single impurity Anderson model. For the particle- 
hole symmetric Anderson model in the KWW approach 
the low energy NRG fixed point is taken to correspond 
to the free conduction chain with the impurity and first 
conduction site removed. This corresponds to V — ¥ oo in 
the bare model, and the equivalent of the strong coupling 
fixed point taken for the Kondo model in the original cal- 
culation of Wilsonii, as the Kondo coupling J ~ V 2 /U. 
The low energy fixed point, however, can equally well be 
viewed as the fixed point of a non-interacting Anderson 
model with a finite hybridization V 0, which we shall 
denote as the Fermi liquid (FL) fixed point. The fixed 
point corresponding to a free NRG conduction chain de- 
pends only on whether the chain has an even or odd num- 
ber of sites. As two sites are removed from the original 
NRG chain in interpreting the fixed point in the KWW 
approach, the KWW and FL fixed points are equivalent. 

When the leading irrelevant corrections to the fixed 
point are taken into account in the KWW calculation, 
the effective Hamiltonian takes the form, 

H N , eff = H NtSC + wxtS N - l),2 5H 1 + 2uj 2 A^ n -^/ 2 SH 2 , 

(39) 
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with 



A = ttV 2 /DA a , where 



5Hi 



E 



(4, 



cr L 2,er 



), «I 2 = £> 1)CT - l) 2 , 
(40) 

where H^ : sc is the Hamiltonian for the truncated free 
conduction chain, 5H\ corresponds to a correction to the 
hopping matrix element between the first and second con- 
duction sites of the truncated conduction chain, and 5H2, 
an on-site interaction at the first site of the truncated 
chain. The corrections to the free energy are then cal- 
culated to first order in uj\ and w 2 , and the following 
expression derived for the impurity specific heat coeffi- 
cient, 



D 



1 + A- 



InA 



(41) 



in units of 27r 2 /ej|/3, and for the susceptibility, 



X = — FT 



a ai uj-2 



D Vl + A-V InA D yl + A -1 / (InA) 2 

(42) 

in units of (^Mb) 2 , where 



a 2 , = 1(1 -A- 1 ), a\ = \(l-k-% 



(43) 



The values of ui\, u>2 in the KWW approach have to 
be calculated from the asymptotic approach of the en- 
ergy levels to their fixed point values. If E P (N) is the 
energy of the lowest single-particle excitation from the 
ground state in the NRG calculation for a chain length N, 
and E* is the corresponding value at the fixed point, the 
value of uj± can be obtained by matching the difference 
E p (N) — E* to that obtained from fixed point Hamilto- 
nian with the leading order correction term SH\ for large 
N. Similarly ui2 can be calculated from the difference be- 
tween the lowest two-particle excitation E pp {N) and the 
lowest two single-particle excitations E P (N) for large N. 

We contrast this with the FL analysis where the two 
leading order corrections correspond to the effective hy- 
bridization with the impurity V and the effective on-site 
interaction U at the impurity site. The hybridization 
term, though a leading irrelevant term, it is technically 
a dangerously irrelevant term, because the fixed point 
changes when V — from that corresponding to an even 
(odd) chain to that for an odd (even) one. The value 
of V can be obtained by requiring E p (N) for large N 
correspond to the lowest single particle excitation of a 
non-interacting Anderson model with hybridization pa- 
rameter V, as described briefly in section IV. The corre- 
sponding value of A for the continuum model is given by 



A A 



-InA 



1 + A- 



1- A" 



(44) 



is a correction factor to the bandwidth D on using 
a discretized NRG chain with discretization parameter 
A(> 1). Once V has been determined U (and J) can be 
calculated from the difference between the lowest two- 
particle and the lowest two single-particle excitations in 
a similar way to the calculation of L)2- Once the renor- 
malized parameters have been determined the physical 
quantities, such as the specific heat coefficient 7, and 
the spin and charge susceptibilities, can be calculated by 
substituting into the relevant RPT formulae. For further 
details we refer to reference o 12 ' 25 . 

We can establish a connection between the two ap- 
proaches using the equations given in Appendix E of the 
KWW I 26 . In this Appendix wi and 0J2 are calculated 
to first order in U for the symmetric Anderson model. 
In the RPT the exact results for 7 and x correspond to 
such a calculation but with renormalized parameters so 
we can use the results given in equations (5.38) and (5.39) 
in KWW I by replacing A (or T in the notation used in 
the paper) by A, and U by U, which gives 



l r l/ 
-2^ 



+ A- 1 )] 



.(InA) 2 D 
alai 2ttA' 



^2 = [-(1+ A 



an 



D 



2ttA 



U_ 
2D 



(45) 



(46) 



We see that both wi and LO2 diverge as A —> 0, if U /irA 
remains finite, which explains why the renormalized pa- 
rameters we calculate on the approach to the QCP tend 
to zero, and the corresponding terms on using the KWW 
approach diverge. 

On using Eqns (|4T|) to (|44)) we find for the continuum 
model, 



2^1 W 2 



3ttA 



2ttA 



U 

7rA 



(47) 



which is the RPT result for the symmetric model. 

There is a clear advantage in the Fermi liquid anal- 
ysis for the asymmetric Anderson model as one deals 
with the same low energy fixed point with just the addi- 
tional renormalized parameter e^, so that results in Eqn. 
(|4T)l are generalized by replacing 1/ttA by p(0). The 
corresponding KWW analysis requires consideration of a 
number of low energy fixed points, strong coupling, inter- 
mediate valence, etc, to cover the full parameter range. 
The generalized results for Eqn. (|47|) . however, can be 
deduced from Appendix D in KWW Il2£ with the sub- 
stitutions, U — > U, r — > A and + U/2 — > id in Eqns 
(D.23-25). 
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